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^ . Abstract 

A new model of dark energy namely "ghost dark energy model" has recently been 
Qi suggested to interpret the positive acceleration of cosmic expansion. The energy 

density of ghost dark energy is proportional to the bubble parameter. In this paper 
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we perform the statefinder diagnostic tool for this model both in flat and non-flat 
universe. We discuss the dependency of the evolutionary trajectories in s — r and 
q — r planes on the interaction parameter between dark matter and dark energy as 
well as the spatial curvature parameter of the universe. Eventually, in the light of 
SNe+BAO+OHD+CMB observational data, we plot the evolutionary trajectories 
in s — r and q — r planes for the best fit values of the cosmological parameters and 
compare the interacting ghost model with other dynamical dark energy models. We 



X: 

■ show that the evolutionary trajectory of ghost dark energy in statefinder diagram 

is similar to holographic dark energy model. It has been shown that the statefinder 
location of ACDM is in good agreement with observation and therefore the dark 
energy models whose current statefinder values are far from the ACDM point can be 
ruled out. 
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I. INTRODUCTION 

Nowadays it is strongly believed that our universe expands under an accelerated expan- 
sion. The various cosmological data gathered from SNe la WMAP SDSS ^ and 
X-ray [4] experiments have provided the main evidences for this cosmic acceleration. Within 
the framework of standard cosmology, a dark energy component with negative pressure is 
responsible for this acceleration. Up to now many theoretical models have been proposed to 
interpret the behavior of dark energy. The first and simple candidate is the Einstein's cosmo- 
logical constant with the time - independent equation of state wa = — 1- The cosmological 
constant suffers from tow deep theoretical problems namely the "fine-tuning" and "cosmic 
coincidence" . In addition to cosmological constant, dynamical dark energy model with time- 
varying equation of state have been investigated to interpret the cosmic acceleration. The 

[9| and dilaton [10] together with interacting dark energy models such as holographic [11 1 
and agegraphic [l2j models are the examples of dynamical dark energy models. The inter- 
acting dark energy models have been constructed within the framework of quantum gravity, 



by introducing the new degree of freedom or by modifying the theory of gravity 13l-ll5|. 
Recently, the Veneziano ghost dark energy has been attracted a deal of attention in the 
dynamical DE category. The Veneziano ghost is proposed to solve the U{1) problem in low- 
energy effective theory of QCD [16| and has no contribution in the fiat Minkowski spacetime. 
In curved spacetime, however, it makes a small energy density proportional to Aqqj^H, where 
^QCD is QCD mass scale and H is Hubble parameter. This small vacuum energy density can 
be considered as a driver engine for evolution of the universe. It is worthwhile to mention 
that this model is totally arisen from standard model and general relativity. Therefore one 
needs not to introduce any new parameter or new degree of freedom and this fact is the 
most advantages of ghost DE. With Aqqd ~ lOOMev and H ~ 10~^^ef, the right order of 
observed DE density can be given by ghost DE. This numerical coincidence also shows that 



this model gets ride of fine tuning probler n |l7l . \l^ Many authors have already suggested 
DE model with energy density as p = aH |l9| . 

Recent observational data gathered from the Abell Cluster A586 support the interaction 
between dark matter and dark energy ^]. However the strength of this interaction is not 



clearly identified 21|. 
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Since many theoretical dark energy models have been proposed to explain the accelerated 
expansion of the universe, therefore the sensitive test which can differentiate between these 
models is required. The Hubble parameter, H = a/a, (first time derivative) and the decel- 
eration parameter q = —aH'^/a (second time derivative) are the geometrical parameters to 
describe the expansion history of the universe, d > or if > means the expansion of the 
universe. Also d > 0, i.e. (? < 0, indicates the accelerated expansion of the universe. Since 
the various dark energy models give ii > 0, g < at the percent time, the Hubble and 
deceleration parameters can not discriminate dark energy models. For this aim we need a 
higher order of time derivative of scale factor. Sahni et al. 22| and Alam et al. [23'], by 
using the third time derivative of scale factor, introduced the statefinder pair {s,r} in order 
to remove the degeneracy of H and q at the present time. The statefinder pair is given by 

'a r — 1 



(1) 



a/73' 3(g-l/2^ 

Depending the statefinder diagnostic tool on the scale factor indicates that the statefinder 
parameters are geometrical. The scale factor a{t) can be expanded near the present time to 
as follows 

ait) = 1 + H^{t - to) - IqoH^it - to)' + ]:roH^{t - to)' + - (2) 

where we consider a(to) = 1 and Hq, qo, tq are the present values of the Hubble parame- 
ter, deceleration parameter and former statefinder parameter, respectively. Up to now, the 
various dark energy models have been studied from the viewpoint of statefinder diagnostic. 
These models have different evolutionary trajectories in {s, r} plane, therefore the statefinder 
tool can discriminate these models. The well known ACDM model is related to the fixed 
point {s=0,r=l} in the s — r plane {22!. The distance of the current value of statefinder pair 
{■So, ''"0} for a given dark energy model from the fixed point {s=0,r=l} is a valuable criterion 
to a model. In addition, the distance of current statefinder values of a given dark energy 
model from the constrained observational value {so,'ro} is a good tool to test a model. 
The dynamical dark energy models that have been investigated by statefinder diagnostic 
tool are: 

the quintessence DE model 



graphic dark energy models 



22 
26 



23 1 , the interacting quintessence models [2J, |25|, the holo 



27j , the holographic dark energy model in non-flat universe 



the phantom model [29|, the tachyon [30|, the generalized chapl ygin ga s model |3l|, the 



interacting new agegraphic DE model in fiat and non-fiat universe 



32 



33|, the agegraphic 



4 



dark energy model with and without interaction in flat and non-flat universe 



new holographic dark energy model 
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J5|, the 
3. 



and the interacting polytropic gas model 
In this work we investigate the interacting ghost dark energy model by statefinder diagnostic 
tool. The statefinder can be applied to diagnose different cases of the model, including dif- 
ferent model parameters and different contributions of spatial curvature. First, we perform 
the statefinder diagnostic in flat universe in sect. II, then we generalize our work to the non 
flat universe in sect. III. In sect. IV, the statefinder diagnostic has been discussed based on 
recent observational data. This work is concluded in sect. V. 



II. INTERACTING GHOST DARK ENERGY MODEL IN FLAT UNIVERSE 

Let us first consider the interacting ghost dark energy in the flat Friedmann-Robertson- 
Walker (FRW) universe. The corresponding Friedmann equation in this case is given by 

H' = ^{Pm+PA) (3) 

where H and ikL are the Hubble parameter and the reduced Planck mass, respectively. 



The energy density of ghost dark energy is given by [50 1 

Pa = aH (4) 
where a is a constant of the model. The dimensionless energy densities are defined as 

o - - o - Pa - P^ 

™ " " 3M2i72' - - 3M^H^ 

Using ([5]), the Friedmann equation ([3]) can be written as 

n„ + = 1. (6) 

In a universe dominated by interacting dark energy and dark matter, the total energy density, 
P = Pm + Pa, satisfies the following conservation equation 

p + 3H{p + p) =0 (7) 



However, by considering the interaction between dark energy and dark matter, the energy 
density of dark energy and dark matter does not conserve separately and the conservation 
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equation for each component is given by 

pm + 'iHpm = Q, (8) 

Pa + 3H{pa+pa) = -Q, (9) 

where Q represents the interaction between dark matter and dark energy. It is worth noting 
that in equation ([8]) the right hand side of ([8]), same as left hand side, should be as a function 
of inverse of time. The simple choice is that the interaction quantity Q can be considered 
as a function of Hubble parameter H such as one of the following forms: (i) Q oc Hp\, (ii) 
Q oc Hpm and (iii) Q oc H{pm + Pa). One can assume the above three forms as Q = Tp^, 
where for case (i) T = 3b'^H, for case (ii) T = ?,lP'H^ and for case (iii) V = Sb'^H-^. 
The parameter 6 is a coupling constant indicating the strength of interaction between dark 
dark energy 5^. The interaction between dark energy and dark matter is also 



matter anc. 
studied in 



53|. Here we assume the third form of interaction for Q. 



Taking the time derivative from both side of Friedmann equation ([3]) and using (El El E]) as 
well as the relation p\ = waPa, one can obtain 

;^ = -^[1 + ^a^a] (10) 

Inserting the third form of interaction term Q = Tp\ = 3b^H-^p\ in the right hand side of 
dH]) and using the relations ffTOj) . the equation of state for interacting ghost dark energy 
in the flat universe can be obtained as 

-1 ^ 2b^ ^ 

In the limiting case of non-interacting fiat universe (i.e., 6 = and VLk = 0), Eq. (fTTl) reduces 
to 



which is in agreement with 5l|. At the early time when « 1, we can see wa = — 1/2 
and at the late time when f^A ~ 1, one can see wa = —1- Therefore the ghost dark energy 
mimics the cosmological constant at the late time. The evolution of EoS parameter of ghost 
model has been studied in [si]]. It has been shown that the interacting ghost dark energy 
model can cross the phantom divide for 6^ > 0.1. 

Using ( ITOl) . the deceleration parameter q in this model can be obtained as 

H 1 3 
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It is clear that at the early time ( when Qj>^ 0) we have q = 1/2 which is equal to the 
value of deceleration parameter obtained in CDM model. Therefore in ghost model, the 
decelerated expansion phase [q > 0) at the early time can be achieved. At the late time ( 
when ~ 1 and wa = —1), we see that g = — 1, which represents the accelerated expansion 
(g < 0) in dark energy dominated universe, as expected. 

Tacking the time derivative of dark energy density parameter in ([5]) and using the ghost 
dark energy density (jl]), we have 



Using ( IT^ and ^Ia = HVt'^ yields 



Q'^ = hlA{l + WKnA) (15) 

where prime denotes the derivative with respect to In a. Tacking the time derivative of ( ITOl) 
and using ([5]), ([9]) and ^ we obtain 

H 9 3 18 

— = -wa^a{wa^a + 3) - -^aw'a + — (16) 

We now find the statefinder parameters {s, r} for the interacting ghost dark energy model 
in the flat universe. From the definition of q and H, the parameter r in ([1]) can be written 

as 

^ = ;^-3g-2. (17) 

Substituting the relations ( fT3l) and (fT6l) in ( fTTl) . the parameter r is obtained as 

9 3 

r = 1 + -WA^AiwA^A + 1) - 2^A^A (1^) 

Inserting Eqs. ( !T3|) and ( ITSjl in the parameter s of ([1]) obtains 

At the late time ( when Qa ~^ 1 ), by inserting wa = — 1 and therefore = 0, the relations 
( |T8l) and ( |T9l) reduce to the constant values (r = 1, s = 0) which refers the statefinder 
parameters of standard ACDM model in the fiat universe. Therefore, from the viewpoint of 
statefinder diagnostic, the ghost dark energy mimics the cosmological constant at the late 
time. 
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By numerical solving of Eqs. ( ITSjl and ( !T9|) . we obtain the evolutionary trajectory of in- 
teracting ghost dark energy in the statefinder plane. It should be noted that in Eqs. f lTSj) 
and ( fT9l) the evolution of i^a and are governed by Eqs. ( ITSl) andf fTTl) . respectively. In 
statefinder plane, the horizontal axis is defined by the parameter s and vertical axis by the 
parameter r. In this diagram, the standard ACDM model corresponds to the fixed point 
{r = l,s = 0}. 

In Fig.(l), we plot the evolutionary trajectories of ghost dark energy model in the fiat uni- 
verse in s — r plane for different illustrative values of interaction parameter b. Here we adopt 
the current values of cosmological parameters Qa and Qm as 0.7 and 0.3, respectively. The 
standard ACDM fixed point {r = 1, s = 0} is indicated by star symbol in this diagram. 
The colored circles on the curves show the present values of statefindr pair {so,ro}. By 
expanding the universe, the trajectories in s — r plane start from right to left. The param- 
eter r decreases, then increases to the constant value r = 1 at the late time. While the 
parameter s deceases from the positive value at the early time to the constant value s = 
at the late time. Different values of interaction parameter b result the different evolutionary 
trajectories in s — r plane. Hence the statefinder analysis can discriminate the interacting 
ghost dark energy model with different interaction parameter. For larger value of b, the 
present values of Sq and tq decreases. The distance of the point {sojTq} form the ACDM 
fixed point {s = 0,r = 1} becomes larger for larger values of interaction parameter b. Fig.(l) 
also shows that the interacting ghost dark energy model mimics the ACDM model at the 



late time. This behavior of ghost dark energy is similar to the holographic new 
agegraphic {32, 33|, chap ly gin gas {stI, generalized chaplygin gas jsil and yang- mils 38 1 



models of dark energy in which they also mimic the ACDM mode^ 



at the late time. 



Unlike the above models, the agegraphic dark energy model 3J, |35|] and polytropic gas 
model 39| mimic the ACDM model at the early stage of the evolution of the universe. The 



evolutionary trajectories of holographic dark energy under granda-Oliveros IR cut-off (new 
holographic model) 36| and also tachyon dark energy model 30|] in s — r plane pass through 
the ACDM fixed point at the middle of the evolution of the universe. The other interesting 
note is that the evolution of ghost dark energy model in s — r plane is similar to the evolution 
of holographic model of dark energy with the model parameter c = 1 in this plane (i.e., see 
Fig. (3) of |27| and upper panel of Fig.(l) in |28| ) . 

Also, it is of interest to discuss the dynamical behavior of ghost dark energy in g — r plane. 
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In g — r plane, we use the geometrical quantity q instead of the parameter s at the horizontal 
axis. In Fig. (2), by solving Eqs. f|T3l) and f|T8l) . the evolutionary trajectories of ghost dark 
energy in flat universe is plotted for different values of interaction parameter 6 in g — r plane. 
Same as statefinder analysis, the q — r analysis can discriminate different dark energy mod- 
els. By expanding the universe, the trajectories start from right to left. The parameter r 
decrease, then increases to the constant value r = 1 at the late time. While the parameter q 
decreases from the positive value ( indicating the decelerated expansion) at the early time to 
the negative value (representing the accelerated expansion) at the late time. Here we see the 
different evolutionary trajectories for different interaction parameters h. The current value 
{go, '"o} can also be affected by interaction parameter. Increasing the interaction parameter 
h causes both the parameters r and g becomes smaller. 




FIG. 1: The evolutionary trajectories in s — r plane for interacting ghost dark energy model in the 
flat universe with the cosmological parameters 17^0 = 0.3 and JIaq = 0.7. The location of standard 
ACDM fixed point is indicated by star symbol. The colored circle points are the location of present 
values of statefinder pair {sq^vq} for different interaction parameter as described in legend. 
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FIG. 2: The evolutionary trajectories m q — r plane for interacting ghost dark energy model in the 
flat universe with the cosmological parameters QmO = 0.3 and Qao = 0.7. The colored circle points 
are the location of present values of statefinder pair {qoj^o} for different interaction parameter as 
described in legend. 



III. INTERACTING GHOST DARK ENERGY MODEL IN A NON FLAT 

UNIVERSE 

In this section we generalize our work in previous section to the non flat universe. The 
Friedmann equation in this case is given by 



10 



where k = 1,0,-1 is a spatial curvature parameter corresponding to the closed, flat and 
open universe, respectively. The dimensionless energy densities of dark energy and dark 
matter have been deflned in ([5]) and dimensionless energy density corresponding to the 
spatial curvature is given as ^Ik = ^^27^- Therefore the Friedmann equation f l2U]) in terms of 
dimensionless parameters is written as 

fi^ + = 1 + ilk. (21) 

Same as previous section, here in the non flat universe, we consider the third form of interac- 
tion between dark matter and dark energy Q oc H{pm + Pa)- Using Eqs. fl20|) and fl2Tl) . this 
form of interaction in non flat universe can be written a.s Q = FpA, where F = 3b'^H^Y~- 
Taking the time derivative of both side of Friedmann equation ( l20i) and using ( |2T1 [HI [9]) as 
well as the relation = waPa, one can obtain 

H 3 

Qk - -[I + + waQa] (22) 



^2 ^2' 



where flk is given by 



1 nl,2 Q -, 

l + ^(H-fi.)(l-=^)-y] (24) 



^4 = — : , 7 = 7^ — (.23) 

Inserting the interaction term Q in the right hand side of continuity equation and using 
the relations @, (!22|) . the equation of state for interacting ghost dark energy in the non flat 
universe can be obtained as 

2 r 1 nh.2 

In the limiting case of flat universe (i.e., Qk = 0), Eq.(!24l) reduces to (fTTl) . as expected. 
Using (1221) . the deceleration parameter q in non flat case can be obtained as 

H 1 3 

g = -1 - = 2(1 + ^'^) + 2^^^^ (25) 

The evolution of dark energy density in non flat universe is obtained by tacking the time 
derivative of (E]) and using the ghost dark energy density (jl]) 

Using ( l25l) and Qa = HQ'j^ results 

Q'j, = Qa{1 + q) (27) 
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where q is defined in (125|) . Tacking the time derivative of Eq. (122|) and using ([5]), ([9]), ( 123|) 
and ([24]) results 

^ 9 3 , ^ ^]A(^]fc + 7)-lo , is 

— = -WA^]A(WA^]A + 3) - -QAiw^ - QkWA _ ^ ) + ^^^^ _ + ^ (28) 

Inserting Eqs. ( 125|) and (!28|) in Eq. (ITTl) . the former statefinder parameter r for interacting 
ghost dark energy in the non flat universe is obtained as 

r = 1 + ^^A^^AK^^A + 1) - ^i^aK - ^fe^A ^n'T ) + ^^^^47?r^Hf^ (29) 

4 z sja — i 4(siA — J-J 



Following 55|, we consider the parameter s in the non flat universe as follows 



(30) 



3(g - nt/2) 

where fi^ = 1 + is a total energy density as defined in Friedmann equation. Obviously, 
in the hmiting case of flat universe, i.e., flk = 0, the above definition is reduced to ([T]). 
Substituting Eqs. ([25]) and ([29]) in ([30]) gets 



,1 + WaQa) - + -TTT^ 7 - 3/2 + — (31) 

2' 3wA S^Qa-I) V GwA J ^ 

In the limiting case of flat universe, the above equations for the statefinder parameter {s, r} 
are reduced to those obtained in previous section. Here in this section, we focus on the 
contribution of spatial curvature on the evolution of ghost dark energy in the s — r and q — r 
planes. For this aim we need to solve numerically the relations ( 125] [29] and [3T]) . Note that 
in these equations the dynamics of EoS parameter wa, density parameter f^A and spatial 
curvature parameter fi^ are given by ([2^ . ([5]) and ([23]) . respectively. 

In Fig. (3), we plot the statefinder diagram for different contribution of spatial curvatures. 
The selected curves are plotted by fixing Q^o = 0.30, Qao = 0.70 and varying Qko = 0.02, 
Qko = 0.00 and Qko = —0.02 corresponding to the closed, flat and open universe, respectively. 
A closed universe with a small positive curvature ( = 0.02) is compatible with some 



observations 



56| . Here we ignore the interaction between dark matter and dark energy and 



focus only on the effect of contribution of spatial curvature on the evolution of trajectories in 
statefinder plane. By expanding the universe, the trajectories evolve from right to left. The 
parameter r decreases, then increases and reaches to the constant value r = 1 at the late time. 
The parameter s decreases forever. The different contributions of spatial curvature exhibit 
the different features in the s — r plane. The colored circles on the curves are the today's 
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value of {scTq} for different spatial curvatures. One can see that the today's value {sojTo} 
of interacting ghost dark energy with different spatial curvatures is discriminated. We can 
clearly identify the distance from a given dark energy model to the standard fiat-ACDM 
model by using the r(s) evolution diagram. Fig. (3) shows that in the closed universe the 
distance of the present value {sq, tq} from the location of ACDM fixed point {s = 0, r = 1} 
is shorter compare with other spatial curvatures. The holographic dark energy model from 
the viewpoint of statefinder diagnostic analysis in the non fiat universe has already been 
investigated in j28|. By comparing Fig. (3) with upper panel of Fig.(l) of 28|, we see the 
similarity of evolutionary trajectories between ghost dark energy model and the holographic 
model of dark energy (with the model parameter c = 1) in non fiat universe. 
Fig. (4) shows the evolutionary trajectories of interacting ghost dark energy in g — r plane for 
different contributions of spatial curvature of the universe. By expanding the universe the 
trajectories evolve from right to left, the parameter r first decreases then increases and the 
parameter q decreases from the positive value at the early time (indicating the decelerated 
phase of expansion) to the negative value at the at late time ( denoting the accelerated 
phase). In g — r plane, the interacting ghost model is discriminated for different contribution 
of spatial curvatures. The current value of statefinder pair {qo,ro} is also distinguished in 
different spatial curvatures of the universe. The value of {go^'^o} is larger in closed universe 
{Qk = 0.02) compare with fiat {Q^ = 0.00) and open {Q^ = —0.02) universe. 



IV. INTERACTING GHOST MODEL AND OBSERVATIONAL CONSTRAINS 

It is clear that constraining the parameterized model against the observational data is 
model dependent. Hence some doubts usually remain on the validity of the constraints on 
the derived quantities such as the present day values of the deceleration parameter and 
the age of the universe. In order to solve this problem, we use the cosmography, i.e. the 
expansion of scale factor in Taylor series with respect to the cosmic time. For this aim, the 
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FIG. 3: The evolutionary trajectories of ghost dark energy model in s — r plane for different 
contributions of spatial curvatures Qko = 0.02 (closed universe), Clko = 0.00 (flat universe), Qko = 
—0.02 (open universe). Here we set ClmO = 0.3 and Oaq = 0.7. The colored circle points are the 
location of present value {sq, tq} for different spatial curvature as indicated in legend. The location 
of ACDM fixed point has been shown by star symbol. 



following functions 



if = (32) 

adt ^ ^ 

adt^ ^ ^ 

I = i^/f- (36) 

adt^ ^ ^ 
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FIG. 4: The evolutionary trajectories of ghost dark energy model in q — r plane for different 
contributions of spatial curvatures ^Iko = 0.02 (closed universe), Qko = 0.00 (flat universe), Qko = 
—0.02 (open universe). Here we set ClmO = 0.3 and Oaq = 0.7. The colored circle points are the 
location of present value {qo,ro} for different spatial curvature as indicated in legend. 



which are namely the Hubble, deceleration, jerk, snap and lerk parameters, respectively 
are introduced. The present values of the above parameters can be used to describe the 
evolution of the universe. For example, qo < indicates the current accelerated expansion of 
the universe and also tq allows to discriminate between different dark energy models. Using 
the Union2 SNela data set jsoj] and the BAO data from the analysis of the SDSS seventh 



release 



60| adding a p rior on h from the recent determination of the Hubble constant by 



the SHOES team |61[ and the age of passively evolving galaxies |58[, the present values 
of the above cosmographic parameters are constrain observationally by using the Markov 
Chain Monte Carlo method 57|]. The best fit values of the cosmographic parameters are: 
{h = 0.718, go = -0.64, tq = 1.02, ko = -0.39, Iq = 4.05} (see table I of [sTI for more 
details) . 
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Inserting the present values of go = —0.64 and tq = 1.02 in Eq. ([T]), the present value of 
statefinder parameter s is obtained as Sq = —0.006. Therefore, observationally, the best fit 
value of the current statefinder pair is {sq = — 0.006, ro = 1.02}. In this section we compare 
the present value of statefinder parameters {s,r} of interacting ghost dark energy that has 
been constrained observationally in 511] with the above best fit value of current statefinder 
pair. 

For this aim we use the best fit constrained values of the cosmological parameters Qmo = 0.35, 
r^Ap = 0.75 and 6^ = 0.08 in the ghost dark energy model that have recently been obtained in 



5l| by using the data of Supernova type la (SNIa) Gold sample, shift parameter of Cosmic 



Microwave Background radiation (CMB) and the Baryonic Acoustic Oscillation (BAO) peak 
from Sloan Digital Sky Survey (SDSS). In Fig. (5) the evolutionary trajectories of interacting 
ghost dark energy in s — r plane (upper panel) and in g — r plane (lower panel) are plotted 
for the above best fit values of cosmological parameters ^Imo, ^ao and 6^. In s — r diagram, 
the evolutionary trajectory starts from {s = 0.86, r = 0.67} at the past time, reaches to 
the {s = 0.08, r = 0.74} at the present time (circle point) and ended at {s = 0,r = 1} 
at the future. The best fit observational value {sq = — 0.006, tq = 1.02} in fiat universe is 
indicated by red-star symbol in this diagram. In g — r diagram, the evolutionary trajectory 
starts from {q = 0.4, r = 1} at the past ( corresponds to the decelerated expansion of the 
universe), reaches to {g = —0.6, r = 0.74} at the present time and ended at {g = —1, r = 1} 
at the late time ( corresponds to the accelerated expansion). The best fit observational value 
{go = — 0.64, tq = 1.02} is also indicated by red-star symbol in this diagram. 
Now we compare the present value {so,ro} of constrained interacting ghost dark energy 
model with other models of dark energy which have been constrained and discussed from the 
viewpoint of statefinder diagnostic. This comparison includes the interacting ghost model, 
holographic, new holographic and generalized chaplygin gas models of dark energy. These 
models have been constrained by astronomical data of SNe+CMB+BAO+OHD experiments 
and also have been discussed in s—r diagram based on the constrained values of cosmological 
and model parameters. This comparison also includes the standard ACDM model as well as 
the best fit observational value {sq = — 0.006, tq = 1.02} in fiat universe. The holographic 
dark energy model with the constrained values (c = 0.84, Q^no = 0.29, Q^o = 0.02, where 
c is the model parameter of holographic dark energy) obtains the today's statefinder pair 



as {sq = — 0.102, ro = 1.357} 281]. The new holographic dark energy model with the 
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constrained values (Qfo/i^ = 0.0228, n^o = 0.2762, Qko = 0.0305, ^ao = 0.6934, a = 0.8824, 
P = 0.5016, where a and /3 are the parameters of model) results the today's statefinder pair 
as {sq = — 0.13,ro = 1.46} [36]. The generalized chaplygin gas dark energy ( GCG model) 
in the fiat universe with the constrained values {As = 0.76, a = 0.033, Qbh"^ = 0.0233, 
Hq = 69.97, where Ag and a are the parameters of the model) gives the today's statefinder 
pair as {sq = —0.007, rg = 1.026} [3l|- Note that the GCG model is constrained in the 
fiat universe, but other models are constrained in general non-fiat universe. Fig. (6) shows 
the location of the present statefinder pair {sqjTq} for the above constrained models as 
indicated in legend. The standard ACDM model and also the best fit observational value 
{so = —0.006, To = 1.02} in flat universe are indicated by black and red star symbols, 
respectively. One can conclude that the ACDM model {s = 0, r = 1} has a shortest distance 
to the best fit observational value {sq = —0.006, tq = 1.02} compare to other dynamical dark 
energy models. Also, the interacting ghost dark energy model has a shorter distance from 
{so = — 0.006, ro = 1.02} compare with the holographic and new holographic dark energy 
models. Among the dynamical dark energy model, the GCG model has a shortest distance 
from the location of observational value in s — r plane. 
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1.1 




FIG. 5: The statefinder diagrams r(s) (upper panel) and r{q) (lower panel) for interacting ghost 
dark energy model. The evolutionary trajectories are plotted in the light of best fit result of SNe 
+ OHD + BAG + GMB, Qao = 0.75, QraO = 0.35 and 6^ = o.08 . The circle points on the 
curves show the todays value {so,ro}, upper panel, and {qo^ro}, lower panel. For comparison, the 
standard AGDM model has been shown by black-star symbol and the constrained observational 
va.lne of isr\.rr\\ a.nd ion.rn]- are indicated bv red-star .svmbol in these dia.ffra.ms 
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FIG. 6: The present value of {so,ro} in the light of best fit result of SNe + OHD + BAO + 
CMB observations for different dark energy model as indicated in legend. The location of standard 
ACDM model and constrained observational value {so,ro} have been shown by black and red star 
symbols, respectively. 

V. CONCLUSION 

Summarizing this work, we investigated the interacting ghost dark energy model in 
statefinder s — r and q — r diagrams. The statefinder analysis can discriminate the in- 
teracting ghost dark energy model for different values of interaction par ameter as well as 
ihe different spatial curvatures of the universe. Like holographic 



26 



28 



32 



33|, chaplygin gas 



new agegraphic 



371], generalized chaplygin gas [31] and yang-mils [38| models of dark 
energy, the ghost dark energy model mimics the standard ACDM model at the late time. 
The evolution of ghost dark energy model in s — r plane is similar to holographic model 
of dark energy with the model parameter c = 1. Different values of interaction parameter 
obtains the different evolutionary trajectories in s — r and q — r planes. The evolutionary 
trajectories r(s) and r(g) for interacting ghost dark energy model in different closed, fiat 
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and open universe has also been investigated. We have shown that different contribution of 
spatial curvatures give the different evolutionary trajectories in s — r and q — r. The spatial 
curvature can also influence the present value of statefinder parameters {sq, ro} and {qo,ro} 
in these planes. Eventually, we performed the statefinder diagnostic for the interacting ghost 
model constrained by observational data. We conclude that the ACDM model {s = 0, r = 1} 
has a shortest distance to the best fit observational value {sq = —0.006, Tq = 1.02} compare 
with other dynamical dark energy models. Therefore the models of dark energy whose cur- 
erent statefinder values locate far from the ACDM point can be ruled out. The interacting 
ghost dark energy model has a shorter distance from {sq = —0.006, ro = 1.02} compare with 
the holographic and new holographic dark energy models. Among the above dynamical dark 
energy models, the GCG model has a shortest distance from the location of observational 
statefinder pair (i.e., {sq — — 0.006, Tq = 1.02}). 
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